Curvature Effects in Gamma Ray Burst Colliding Shells 
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ABSTRACT 

An elementary kinematic model for emission produced by relativistic spherical colliding shells 
is studied. The case of a uniform blast-wave shell with jet opening angle 9j ^ l/F is considered, 
where F is the Lorentz factor of the emitting shell. The shell, with comoving width Ar', is assumed 
to be illuminated for a comoving time At' and to radiate a broken power-law vL^, spectrum 
peaking at comoving photon energy e^^. g . Synthetic GRB pulses are calculated, and the relation 
between energy flux and internal comoving energy density is quantified. Curvature effects dictate 
that the measured vF^, flux at the measured peak photon energy Cpk is proportional to e^^, in the 
declining phase of a GRB pulse. Possible reasons for discrepancy with observations are discussed, 
including adiabatic and radiative cooling processes that extend the decay timescale, a nonuniform 
jet, or the formation of pulses by external shock processes. A prediction of a correlation between 
prompt emission properties and times of the optical afterglow beaming breaks is made for a 
coohng model, which can be tested with Swift. 

Subject headings: gamma-rays; bursts — gamma-rays: theory — radiation processes: nonthermal 



1. Introduction 

In the coUapsar scenario for GRBs, pulses in 
GRB light curves are thought to be produced by 
collisions between relativistic shells ejected from 
a central engine (see Zhang & Meszaros (2004) 
for a recent review). The interception of a more 
slowly moving shell by a second shell that is 
ejected at a later time, but with faster speed and 
larger Lorentz factor, produces a shock that dis- 
sipates internal energy to energize the particles 
that emit the GRB radiation. This scenario is 
widely considered to explain pulses in GRB light 
curves (Kobayashi, Piran, & Sari 1997; Daigne & 
Mochkovitch 1998). Studies of pulses are impor- 
tant to decide if GRB sources require engines that 
are long-lasting or impulsive (Dermer & Mitman 
2003), with important implications for the nature 
of the central engine, which is often argued to be a 
newly formed black hole powered by the accretion 
of a massive dense torus. 

Here we construct an elementary kinematic 
model for colliding shells, assumed spherical and 



uniform within jet opening angle 9j. This is the 
sort of jet that Frail et al. (2001) discuss regard- 
ing the standard energy reservoir result, where jet 
opening angles are inferred from the time of achro- 
matic spectral breaks in optical afterglow light 
curves. 

We also perform this study in order to quantify 
the curvature constraint of a spherically emitting 
shell traveling with bulk Lorentz factor F, which 
implies that the shell radius 



2T'^cL, 



rl{l + z) 



(1) 



in order to produce variability on timescale t^jar 
(Rybicki & Lightman 1979; Fenimore, Madras, & 
Nayakshin 1996). This study also quantifies the 
rate at which fiux decays at a given energy due to 
curvature effects, and the range of validity of the 
approximate relation 



cr^u'^V^/dl 



(2) 



between internal comoving energy density Uq and 
observed energy fiux where is the luminos- 
ity distance (See Appendix A). The accuracy of 
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this relation is important to quantify 77 opacity 
constraints (Lithwick & Sari 2001; Dermer 2004) 
applied to GRB pulses as measured with the GRB 
monitor and Large Area Detector on GLAST^, 
as well as to make estimates of photomeson pro- 
duction in GRB blast waves (Waxman & Bahcall 
1997). 

If curvature effects dominate the late-time emis- 
sion in GRB pulses, then a unique relation is found 

whereby the value of the i^F^ peak flux f^^^, (in cgs 
units of ergs cm~^ s~^) at peak photon energy Cpk 
decays in proportion to oc e^^. This relation is gen- 
erally not observed in long, smooth GRB pulses 
studied by Borgonovo & Ryde (2001), who find 
power-law decays f^^^ oc e^^,, with 0.6 < C ^ 3. 
Remarkably, values of ( for different pulses within 
the same GRB arc confined to a rather narrow 
band of values. The wide range of values of are 
found not only in multi-peaked GRBs, but also 
in single-peaked GRBs that display smooth fast- 
rise, slow-decay light curves (Borgonovo & Ryde 
2001; Ryde & Petrosian 2002). The smooth sin- 
gle peak GRBs could arise from curvature effects 
(Fenimore, Madras, & Nayakshin 1996), or to ex- 
ternal shocks (Dermer, Bottcher, & Ghiang 1999). 
For GRB pulses that could be produced by spher- 
ically symmetric shell collisions, discrepancy with 
observations suggest a breakdown of our assump- 
tions. 

In the next section, the kinematic model is pre- 
sented. Calculations based on this model are pre- 
sented in Section 3. In Section 4, we discuss the 
possibility that radiative-cooling effects produce 
the power-law relation, implying a prediction that 
can be tested with Swift^. Alternately, the uni- 
form spherical shell assumption could break down, 
or the basic model of colliding shells could be in 
error. The Appendices give derivations of simple, 
widely-used approximations related to this study, 
a derivation of the curvature relation f^^^^ oc e^j,, 
as well as an analytic form for the time-dependent 
pulse profile, leading to a simple expression for 
the light curve of a pulse in the curvature limit. A 
brief summary is given in Section 5. 

^ glast .gsfc.nasa.gov 
swift .gsfc . nasa. gov 



2. Kinematic Model 

A simple kinematic model for the received flux 
from the illumination of a spherically symmetric 
shell resulting from shell collisions is studied. A 
shell with finite width is assumed to be uniformly 
illuminated throughout its volume for a fixed du- 
ration during which the shell travels with constant 
speed from the explosion center. Light-travel time 
and Doppler effects are treated without regard to 
details of the energization and cooling of the ra- 
diating particles. This approach gives kinematic 
expectations of curvature effects in a GRB collid- 
ing shell system. 

The uFi, flux measured at dimensionless photon 
energy e = hu/meC^ and time t is given by 

-1 /'27r />! /.oo 

xe'f{e',^i',cf>':rj') , (3) 

where primes refer to comoving quantities, the in- 
tegration is over volume in the stationary (explo- 
sion) frame, the Doppler factor 



(3 = s/1 - l/r2, and f' = (1 + z)e/S (see Granot, 
Piran, & Sari (1999), noting the correction of a 
(1 -|- z) factor in the relation between the emitted 
and received photon frequencies). The emissiv- 
ity j;(e*,f7*) = dE^/dV^dUdn^de^ = 6^f{e',n'), 
where = SI* is the directional vector (/i, 0), 
fj,' = (/i — /3)/(l — /3/i), and c6' = (p. We use a nota- 
tion where asterisks refer to quantities in the sta- 
tionary frame (though we have dropped asterisks 
for the spatial variables r and f2), and unscripted 
quantities refer to the observer frame. 

The blast wave is assumed to emit isotropically 
in the comoving frame, which could apply to syn- 
chrotron and synchrotron self-Compton processes 
with randomly-ordered magnetic fields and elec- 
tron pitch-angle distributions, but not to external 
Compton processes. Moreover, the observer is as- 
sumed to be located along the azimuthal symme- 
try axis of the jet, or is viewing a uniform jet with 
opening angle 6j 3> l/P. Therefore 

m = ^ f df, p dr rh'j'ie'; r, t') , 

(5) 
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noting that S{r) = S for a uniform jot. The 
emissivity is related to the internal energy density 
Ue'{r,t') through the relation 



e'/(e';r,iO 



Ar' 



(6) 



where Ar' = TAr is the proper shell width, and 
the mean escape time of photons from the shell 
volume is approximated by Ar' /c. 

Further consider a uniform jet with no angular 
dependence other than that the emission goes to 
zero a,t 9 > 0j = arccos/Xj. The emitting shell is 
assumed to be illuminated for the comoving dura- 
tion tg < t' < tg + At'. The spectrum is appprox- 
imated by a broken power law with peak pLi, flux 
at energy e^^, given by the expression 

u,,{r,t')=u'aH{t';t'o,t'o + At') 

x[x''H{l-x) + x''H{x-l)] , (7) 

where H are the Heaviside functions, a(> 0) and 
6(< 0) are the i^Lj, indices, and x = e'/e^j, q = 
{l+z)e/6e'pf. Q. The peak vF^, comoving photon en- 
ergy e'pi^ Q is also supposed to be constant through- 
out the shell. The total integrated photon energy 
density for this spectrum is u^^j = u'Q{a~^ —b~^). 

The observing time t is related to the emitting 
time measured in the stationary explosion frame 
through the relation 



t ru 
f — — + lL 

1-1-2: C 



(8) 



The zero of time = corresponds to the moment 
of shell ejection, with the location of the inner edge 
of the shell given by the relation ri(i,) = fict^. 
The first moment of shell illumination takes place 
when the inner edge of the shell is at radius Tq = 
/3ct*o = /3rctQ for a shell moving with constant 
speed r. 

The finite shell width and finite duration of the 
illumination implies two constraints on the inte- 
grations over r and /i. The shell-width constraint 
j3ctt, <r< (3ct* + Ar implies 



<r < 



Pct^ + Ar 



(9) 



Due to light-travel time effects of the relativis- 
tically moving shell, values of r contributing to 
the signal observed at time extend over a range 



Ar/(1 - |3^l) ~ r^Ar (Rees 1966; Granot, Piran, 
& Sari 1999). 

The illumination constraint t'^ < t' = t*/r < 
I'q + At' implies 

-{'^-ct,)<r<-C-^-ct,+cTAt'). (10) 

The zero of observer time is when a hypotheti- 
cal photon ejected at = and r = from 
the inner edge of the shell would reach the ob- 
server. The time at which the signal is first 
detected by the observer is therefore given by 
^^nit ^ _ _ Ar]/f3c ro/2r2c - (Ar/c), 
and the observing time when a photon emitted 
from the inner edge of the shell at the first in- 
stant of shell illumination reaches the observer is 
t^o = ro(l - /3)//3c ro/2T'^c. The final expres- 
sions in these last two relations hold in the limit 
r » 1. Hence 



x[x''H{l-x)+x''H{x-l)] . (11) 

where r; = inax.[Pctz/{i- - /3n),{ro//3 - ctz)/lj] 
and r„ = rmn[{l3ctz + Ar)/(1 - Pn), (ro//3 - ct^ + 
cTAt')/t,]. 

3. Calculations 

We examine the accuracy of the approximate 
expressions relating r and tyar, eq. (1), and en- 
ergy flux and internal energy density, eq. (2) (see 
Appendix A) . Let t j represent a fiducial variabil- 
ity time scale for the observer. We introduce ra- 
dius, time, and width parameters, denoted by rjr, 
rit, and r]A, respectively, to relate tf to source- 
frame quantities. The curvature constraint for 
blast-wave radius suggests that we write 

r = 2r]rT'^ctf/{l + z) . (12) 

Because dt' = 6dt/{l + z), we define 

At' = 2r7]ttf/{l + z) (13) 

to relate the intrinsic variability time scale At' in 
the comoving frame to tf. The comoving width 
of the emitting region Ar' < At'/c, by causality 
requirements (if it were larger, then large ampli- 
tude variability would not be possible except for 
random statistical fluctuations). Thus we define 



Ar' = 2Tr]Actf/{l + z) , 



(14) 
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Fig. 1. — Light curves of a causal pulse at differ- 
ent dimensionless observing energies, for a model 
with standard parameters (see text) and rjr = r]t = 
r/A = 1- Also indicated by the arrows are light 
curves at e = 0.1, 1.0, and 10 for a model with 
% = VA = ^ and rjr = 0.1. 

with the causal requirement rjA < When 
Ar' <C cAf/ , the duration of the emitting region is 
not determined by its causal size scale, but rather 
by the duration of emission radiated from a region 
much smaller than At' /c. 

We solve eq. (11) for the following standard pa- 
rameters: r = 300; z = 1 (so that dL = 2.02 x lO^^ 
cm for a ACDM cosmology with = 0.27, 
f^A = 0.73, and Hubble's constant of 72 km s~^ 
Mpc~^, as implied by the WMAP results; Spergel 
et al. (2003)); e^^ g = {l+z)epk,o/'2T with epk,o = 1 
(i.e., peak photon energy at the beginning of the 
pulse equal to 511 keV); u'q = 1 erg cm~^, a = 4/3; 
b = —1/2; and tf = 1 second. 

Fig. 1 shows the appearance of a kinematic 
pulse with rjr = % = r]A = i 8± & number of 
photon energies. Also shown in Fig. 1 are kine- 
matic pulses formed when 774 = r/A = 1, and 
rjr = 0.1. Note the characteristic rounded, weakly 
asymmetrical (on a linear scale) light curve shapes 
that are formed when Ar' = cAt'. Time de- 
lays from different parts of the width of the emit- 
ting shell are important to determine the pulse 
shape in this case. The smaller emitting vol- 
ume when the shell is energized at O.lro rather 
than at vq produces a pulse with a fluence smaller 
by a factor f^'-o+^rAt' ^^^y rro+cFAt' ^^^2 ^ 

mT'^" '^^^^ I Iro" ~ (1-1" - 0.l3)/(8 - 1) ^ 

1/7 smaller. Indeed, 1/7 is the asymptotic limit of 



Fig. 2. — Evolution of the spectral energy distri- 
bution due to curvature effects for a model with 
standard parameters (see text). In the declining 
phase of the pulse, the value of fe^,, oc e^j., as in- 
dicated by the dashed line. 



the fluence reductions due to the different volumes 
illuminated for a flare lasting for the equal proper 
times, but in one case emerging from deep within 
the jet, and in the other case where the illumi- 
nation begins at the location ro = 2r^ct//(l -|- z). 
The intrinsic duration of the pulse, when combined 
with the curvature effects, has produced a pulse 
with FWHM duration of » 2 seconds, as compared 
with the fiducial time scale of 1 second. Thus the 
combined width, duration, and (off-axis) curva- 
ture effects have lengthened the basic timescale by 
a factor of about 2 at e = epk,oj with a narrower 
FWHM duration when e > €pk,o and a broader 
FWHM duration when e < epfe,o- 

Fig. 2 shows the evolution of the spectral energy 
distribution for this pulse. Notice the rapid decay 
oc Cpf. of the I'Fi, peak flux /e^j^ measured at epk 
during the decay portion of the pulse. This behav- 
ior is characteristic of all pulses where curvature 
effects from off-axis emitting regions dominate the 
late-time behavior of the light curve. 

Fig. 3 shows characteristic light curves when 
Ar' <C cAt', that is, when the shell is very thin 
compared with the size scale associated with the 
intrinsic pulse duration. These light curves exhibit 
pulses that have much sharper peaks than in the 
general case of Fig. 1, and that are asymmetri- 
cal with a distinct trailing edge of emission. The 
fluence contained in the pulse is the same as in 
the pulse of Fig. 1 with r?A = 1, but the FWHM 
pulse duration at e = 1 is w 1 s, comparable to 
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Fig. 3. — Light curves of a thin-shell pulse at dif- 
ferent dimensionless observing photon energies for 
a model with rjr = Vt = ^ and r/A = 0.1, so that 
the shell width Ar' = O.lAt'/c. 



tf, while the peak flux is about twice as large, 
due to the different geometry. The difference in 
geometries of a causal and thin-shell pulses intro- 
duces a physical effect required for accurate calcu- 
lations of scattering or opacity processes in GRB 
blast waves. When Ar' w cAt', then the pho- 
ton field can be considered to be roughly isotropic 
for scattering and opacity calculations. But when 
the shell radiates for a much longer time than the 
light-crossing time through the width of the shell, 
that is, when At' » Ar'/c, the geometry of the 
outflowing photon flux is much more anisotropic, 
giving higher thresholds and lower rates for 77 and 
photohadronic processes due to the reduction in 
the frequency of head-on collisions. 

Fig. 4 shows the characteristic light curve 
shapes formed when curvature effects dominate 
the temporal evolution of the light curve (Fen- 
imore, Madras, & Nayakshin 1996). Here the 
pulses are very asymmetrical, with a sharp lead- 
ing edge. In this calculation, rj^ — 1 and 
Va = Vt = O-l- The peak i^F^ flux at e = 1 
reaches a value of only 1.4 x 10^^ ergs s~^ with a 
duration of « 0.4 s. The total energy released is 
smaller by a factor of 10 than in Fig. 1 with tja 
and Fig. 2 as a result of the shorter intrinsic pulse 
duration. 

The bottom panel in Fig. 4 shows, in a log-log 
relation, that the flux decays as at e > epk- 

When e < Cpk, the flux decays as at early 

times, breaking to a behavior at late times 

due to curvature effects. 
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Fig. 4. — Light curves of a curvature pulse at dif- 
ferent dimensionless observing photon energies for 
a model with rjr = 1, but with tja = ijt = 0.1. 

The spectral and temporal behavior of the cur- 
vature pulse can be derived in the (5-function ap- 
proximation (not to be confused with the Doppler 
factor S). Letting e'i'(e'; r, t') oc e'"(5(r' - ro)5(t' - 

t'o) = e"'S{r-ro)S[t-{l + z)t'JS] in eq. (5), using 
the invariance of the 4-volume, then 

m oc ^ r° J d{l - /Jm)5^-«<5[(1 - - ^] 



oce°( 



2-a 



(15) 



To ro 
where = ro/fiTc, and a is the energy index. 
(This result corrects the expressions given by Fen- 
imore, Madras, & Nayakshin (1996) and Ryde & 
Petrosian (2002), where the delta- function pulse 
in time, 5{t' — Iq), is not transformed between the 
comoving and observer frames.) The dependence 
in eq. (15) is derived more carefully in Appendix 
B, and analytic forms for the pulse proflle, includ- 
ing a simple functional form for the pulse profile 
in the curvature limit, are derived in Appendix C. 
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Discussion 



spectral energy density is thus 



The estimate L = Andl^E - 4:TrrQCUQT'^ , where 
the received flux is intensified by two powers of T 
for the relativistic time contraction and photon en- 
ergy enhancement in a blast wave geometry, is gen- 
erally used to relate bolometric energy flux and in- 
ternal energy density (Appendix A) . More remark- 
ably, the allowed radius of the radiating spherical 
shell is « 2r^ times larger than inferred through 
causality arguments applied to the measured vari- 
ability time scale. This effect greatly dilutes the 
comoving photon density compared with a station- 
ary emitting region, and essentially explains the 
unusual properties of GRBs. In total, we see that 



(16) 

For the nominal parameters used in the figures, 

$B ^ 4.8 X 10-i42„,(s) u'oT%o ergs cm'^ s-K 

The most rigorous limits on 77 attenuation are 
obtained by determining the minimum value of 
the product UgAt' that can produce a pulse with 
a measured peak flux f^^,^ and full-width half- 
maximum duration ti/2 for a given value of F. It 
is the product u^At' that enters into the 77 atten- 
uation (and photomeson) calculations. If the shell 
is found to be optically thick at some photon en- 
ergy for a given value of F even in this case, then F 
must be larger if photons with the corresponding 
energies are detected. 

Inspection of the various cases shows that the 
pulse formed in the curvature limit produces the 
brightest measured flux and shortest duration for 
a given value of the product u'^At'. This is because 
the measured duration is due entirely to curvature 
effects, and the radiated energy is compressed into 
the shortest duration and brightest pulse in this 
limit. From the results of Appendices C and D, 
this implies that the expression 



u'oAt' ^ 



[2 



l/(3-a) 



8c3F5ti/2 



(17) 



gives the smallest possible values for UqcAI', and 
this expression will therefore yield the most re- 
liable minimum Lorentz factors for 77 attenua- 
tion calculations derived from GLAST or ground- 
based air Cherenkov telescope observations. The 
corresponding expression for the comoving photon 



0.26(1 + ^)4/,^, 
8c^rHi/2At' 



(18) 

where ti/2 is determined at photon energies near 
the peak of the uFi, spectrum. Eq. (18) is a factor 
of 3 smaller than the expression used for a comov- 
ing spherical blob with 5 — > F and tvar — * ti/2 (see 
Eq. [2] in Dermcr (2004)). 

Three generic types of pulses have been iden- 
tified for the simple kinematic pulse, namely the 
curvature case where ro ^ cTAt' , the causal case 
where ro ~ FAr' « cTAt' , and the thin shell case 
where Ar' <C cAt'. In all three types of kinematic 
pulses, curvature effects dominate the formation 
of the spectrum at late time t J> {l + z)At'/2T, so 
that /epj. oc Epj. if curvature effects dominate pulse 
formation at late times. 

The curvature relationship can be derived from 
a simple scaling argument by noting that the 
differential stationary-frame shell volume which 
contributes to the received flux, given by dV = 
2TTrQArdfi, remains constant with time. This is 
becaiisc the relation between reception time t and 
/U for a shell that is instantaneously illuminated 
at comoving time is t = (1 -I- z)TtQ{l — j3^i), so 
that djjL oc dt. The vF„ flux /, = S^L'/A-ndl = 
6'^V'e'j{e')/4Trdl = 6^Ve'j{e')/ATrdl, where L' is 
the comoving luminosity of the emitting volume 
that contributes to the flux at time t. For an emis- 
sion spectrum that is flat, that is, e'j(e') oc e'°, 
/epfc <^ ^pk because epk oc i5 for a uniform shell. 

Analysis of BATSE GRB light curves (Bor- 
gonovo and Ryde 2001) shows that the peak fltixes 
of a GRB pulse generally follow a relation whereby 



pk ^ ^pk 



(19) 



Values of ( for different GRBs vary over a wide 
range from w 0.6 to 3, with values of C roughly 
constant for piflscs within the same GRB or in a 
GRB consisting of a single smooth pulse. In most 
GRBs, therefore, curvature effects do not make a 
large contribution to the decay phase of a GRB 
light curve. 

An interesting question is the source of the dif- 
ference of observations from our kinematic model 
pulses. One possibility is that the jet has angular 
structure, and varies with directional energy re- 
lease and baryon- loading on angles 6 w fewxF"^. 
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The anglc-dcpcndcnt speeds in such a system 
would produce a deformed coUiding shocked fluid 
shell where the spherical symmetry assumption 
fails, as therefore would the uniform jet model. 
If this is the case, then GRB prompt emission 
data can in principle be analyzed to reveal shell 
structure and to determine whether this behavior 
is consistent with a universal jet structure (Zhang, 
Dai, Lloyd- Ronning, & Meszaros (2004); see Frail 
(2003) for review). 

Rather than treat these geometrical effects 
here, we consider instead whether radiation effects 
could form a power-law relation between /e^^ and 
Cpfc. The most naive system considers a fixed vol- 
ume of shocked fluid within which the peak of the 
uF,^ spectrum is made by a large population of 
quasi-monoenergetic electrons that radiates most 
of the power through the synchrotron process in 
a mean magnetic field of strength B. If these 
electrons mainly have comoving Lorentz factors 7, 
then their luminous power oc S^7^. Because the 
peak of the i^F,y spectrum is oc Bj"^, and assuming 
B is constant, then f^^^ oc Cpk or (gyn = 1 for this 
simple synchrotron model with constant magnetic 
field. This model can therefore only apply in rare 
cases. 

A better treatment must consider the evolution 
of 7 due to synchrotron and adiabatic losses in the 
expanding shell. The equation of electron energy 
evolution is given by 



Eq. (20) becomes 



"dtf 



1 



dt' 3 Q-KTUeC 



(20) 



where the comoving shell volume changes with 
time according to V^'^ oc t'^™, with m = corre- 
sponding to no expansion, and m = 1 correspond- 
ing to 3-dimensional expansion. 

The magnetic field will also change as a re- 
sult of the expansion of the shell volume. In 
the flux freezing limit where the magnetic field 
is randomly oriented, BB? oc const, implying 
B oc V'^^^ oc i''2m r^YiQ well-ordered magnetic 
field required to explain the polarization observa- 
tion of the GRB 021206 observed with RHESSI 
(Coburn & Boggs 2003) suggests that there is not 
an efficient mixing and randomization of the mag- 
netic field directions. For simplicity, we therefore 
write B oc t'~'^^™', where v = 1 gives the flux- 
freezing limit. 



dr T 



(21) 



where r > 1 is a dimensionless time variable, and 
vq is a dimensionless synchrotron energy loss rate. 
Eq. (21) is analytic, but it is sufficient to consider 
two limiting cases of dominant adiabatic losses or 
dominant synchrotron losses at late times. In the 
case of dominant adiabatic losses we have (drop- 
ping the primes) 7 oc B oc and 
€pk oc 4-2^(1+-), so that fe.Jepk oc B oc e;{'^+'^ 
Thus Cadi = l + [v/{v + l)]. Evcn for a wide range 
of values oi v, 1 < Cadi ^ 2, and (adi is indepen- 
dent of the geometry factor m. 

If synchrotron losses dominate the cooling, 
-d'y/dt oc 5^72. The dependence B oc t'^''"' 
therefore implies 7 oc t^um-i^ gg ^j^g^^ ^ 



Hence fe^^/epk oc B oc e. 



vm/ {l—3vm) 
pk 



SO 

that (syn = 1 + [vm/{l — 3vm)]. Except when 

V <^ 1, (syn ~ 1 when TO = and Csyn ~ 0.5 - 

0.67 when m = 1. In this simple model, therefore, 
values of C between 1/2 and 2/3 are only possible 
for 3-dimensional expansion. 

Three-dimensional expansion is more likely to 
occur for narrowly coUimated blast waves than for 
blast waves with large opening angles, and the nar- 
rowly collimated blast waves would have "beam- 
ing breaks" in the optical afterglow light curves at 
earlier times. If our conjecture that the ftj,kl^pk 
relationships are due to synchrotron and adiabatic 
effects in GRB blast waves with different opening 
angles, then those blast waves with C < 1 should 
be correlated with earlier beaming break times. 
Because this effect is only seen when synchrotron 
losses dominate the cooling, GRBs with C, < 1 
should also display cooling spectra with photon 
indices « 3/2 below €pk. 

Borgonovo & Ryde (2001) find several GRBs 
and many pulses in the BATSE sample with statis- 
tically significant values of Q less than unity. These 
GRBs however preceded the afterglow era. For 
those GRBs which have measured beaming breaks 
(sec Table 1 in Bloom, Frail, & Kulkarni (2003)), 
only GRB 990123 has sufficiently bright BATSE 
data to provide a data point for such a correlation. 
GRB 990123 has not yet been analyzed to give C, 
while analysis of Beppo-SAX data is in progress 
(F. Ryde, private communication, 2004). 

Such a model for the fe^kl^pk relationship 
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would explain why C is approximately constant 
for different pulses within a GRB, provided that 
the opening angle of the GRB jet remains the same 
throughout the period of activity of the GRB en- 
gine. 

The adiabatic/synchrotron model would not, 
however, explain pulses with 2 < C < 3. There are 
many such pulses in the Borgonovo & Ryde (2001) 
sample, though generally with large error bars. If 
analysis of Beppo-SAX or Swift data reveal such 
GRBs, then another explanation is required. One 
possibility is that GRB pulses are due to the in- 
teractions of a single impulsive blast wave with 
inhomogeneities in the surrounding medium. This 
version of the external shock model for the prompt 
phase can be much more efficient than an inter- 
nal shell model (Dermer & Mitman 1999, 2003), 
and permits quantitative studies of the statistics 
of BATSE GRBs (Bottcher & Dermer (2000); see 
Zhang & Meszaros (2004) for a review of the in- 
ternal/external controversy). 

Predictions for the Upkhvk relationship in 
an external shock model (Dermer, Chiang, & 
Bottcher 1999) can be derived by adapting the 
equations for blast wave deceleration in a uniform 
medium with F = ro/[l + [xlxd)^\, where Fq is 
the initial Lorentz factor, Xd is the deceleration 
distance and g is the radiative index {g = 3/2 
and 3 for an adiabatic and fully radiative blast 
wave, respectively). In the deceleration phase, 
X cx ti/(2s+i) and therefore V oc i-9/(2i?+i). in this 
model, epfe cx ^B^"^^^ and oc V"^ B"^ where 
oc oc t-4g/(2g+i) -^^ slow-cooliug regime, 
and 7pfc oc (arF)"^ oc f-2/(2s+i) in the fast-cooling 
regime. Thus U^^^jtyk oc BT. 

In the slow-cooling regime, epk oc F* oc 
f-4g/{2g+i)^ and /^^^ OC e^f . In the fast-coohng 
regime, Cpk oc f-2/(2s-M) and f,^, oc e^+». In the 
slow-cooling and fast-cooling regimes, therefore, 
values of (sc = 3/2 and C/c = 1+5, respectively, 
are predicted. Provided that the surrounding 
medium is uniform (which can be inferred from 
afterglow modeling, though at a larger distance 
scale), the slow-cooling result implies a definite 
value of (sc = 3/2 for fast-rise, smooth decay light 
curves when spectral analysis demonstrates that 
the GRB evolves in the slow-cooling regime. For 
GRBs in fast-cooling regime, this estimate implies 
5/2 < (fc < 4, and in these cases cooling spectra 



should be apparent. Further work will be needed 
to extend the results to radial density gradients 
of the circumburst medium, and to verify that 
these relations hold for deceleration in small den- 
sity inhomogeneities that form GRB pulses in the 
external shock model. 

5. Summary 

A simple kinematic model for GRB colliding 
shells has been constructed that provides a frame- 
work for analyzing radiative processes in a simpli- 
fied geometry of a thin or thick shell traveling at 
relativistic speeds. The relationship between ob- 
served flux and comoving photon energy density 
for a given value of F has been studied, showing 
that the curvature limit yields the smallest value 
of the product UfyAt'. This result can then be 
used to deduce conservative lower limits on bulk 
Lorentz factors derived from the condition of 77 
transparency. 

The kinematic model predicts the ciirvature re- 
lationship /g^j. oc e^j,, with ^ = 3, at late times in 
GRB pulses. Equivalently, curvature effects im- 
ply that fe{t) oc and epk oc t''^. BATSE 
data for GRB pulses do not display the curva- 
ture relationship in most cases (Borgonovo & Ryde 
2001), suggesting that the physics of pulse for- 
mation is dominated by other effects. A simple 
model for joint evolution of /e^^^ and Cpk that takes 
into account adiabatic and synchrotron losses im- 
plies that 1/2 < C < 2, and that C ~ 0.5 only 
when the shell undergoes three dimensional expan- 
sion and the electrons which produce the emission 
near epk are rapidly cooling through synchrotron 
losses. Spectral analysis of Swift data, and corre- 
lations of C with times of the beaming breaks in 
optical afterglow light curves, can test this predic- 
tion. Such a correlation would validate an adia- 
batic/synchrotron model for GRB prompt radia- 
tion, relate properties of the prompt phase with 
the afterglow, and provide key insights into the 
properties of GRB jets. 

Another possibility is that the /e^j, oc e^^ rela- 
tionship is formed by external shock processes, and 
a simple derivation of C was given for blast-wave 
deceleration in a uniform surrounding medium. 
Analysis of prompt data has the potential to test 
the external shock model, though complications 
regarding density gradients and inhomogeneities 
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in tlic' c ircumburst medium must be considered in 
more detail. 

A final possibility is that jet structure produces 
the measured relationship between f^^^, and epk- 
The validity of a universal jet model will be tested 
by determining whether observed values of ( can 
derive from the proposed angle dependence. In the 
meantime, comparing the predictions of the adia- 
batic/synchrotron and external shock models with 
time-resolved spectroscopy of GRB pulses and af- 
terglows has the potential to rule out or validate 
these models. 

I thank Markus Bottcher, Felix Rydc, and the 
anonymous referee for valuable comments. This 
work is supported by the Office of Naval Research 
and NASA GLAST Science Investigation grant 
DPR S-13756G. 



A. Relations and Estimates 



First wo relate the total energy in photons between the eomoving and stationary frames. The differential 
number of photons iV*(e*,ri*) per unit energy and solid angle transforms as Af^*(e*,ri*) = SN'{e' ,^'), as is 
easily seen by calculating the Jacobian of the transformation, or by noting that e~^dN/dedfl = e~^A^(e, f)) is 
an invariant. For an isotropic, monochromatic photon spectrum in the eomoving frame, iV'(e', fl') — N()6{e' — 
eo)/47r, and the total photon energy in the eomoving frame is just Eq = NqCq (in units of the electron rest 
mass). The differential photon spectrum in the stationary frame is therefore A^*(e*, fi*) = 5Aro5(e*/(5 — eo)/47r 
= d'^No6{e^ - (5eo)/47r, so that 

dn, J de, e* N4e,,n,) = J d/i = TE' . (Al) 
This result is obvious by noting the symmetry of the transformation equation e* = re'(l + fSfj,') with respect 

to /i'. 

Because — L/Andj^, by definition of the luminosity distance dz,, the fluence (p = ^sit) = L^{t)/AiTd\ = 
i*(t*)(l + z)/AiTd\, where (t) and (t*) are times of reception and emission in the observer and stationary 
frame, respectively. Using eq. (Al) gives 

^=^(1 + ^). (A2) 

A simple estimate relating eomoving energy density u'q with energy flux <!>£; is obtained by noting that 
the stationary frame luminosity of a blast wave is given by L* = d£'*/dt* = F^L', where L' = dE'/dt' = 
Uo47rr^Ar'/(Ar'/c). Thus — cr'^u'^T'^ /d\, giving eq. (2). If the variability is produced by curvature 
effects according to eq. (1), then 

{l + zfdl • ^ ' 

Note that the same basic dependence, though with F replaced by 5, is derived for a (eomoving) spherical 
blob geometry. In this case, $b = S^L' /Awd!^, and L' = 47rr^^C'Uo/3, with blob radius = c6tyar/{^ + z). 



B. Analytic Derivation of fupki^) vs. Cpfc Relation 
Starting with eq. (5), we approximate 

e'j'(e'; r, t') = Ke"^6{r' - ro)5{t' - QH{e'; ej, . (Bl) 

Normalizing to the eomoving energy E'^ of a pulse implies that 

_ aE'p 

2^{l-H^)rl{e'S-e'n- ^^^^ 

The integrals can now be performed. First note the subtlety that \dr/dr'\ = 6, whereas Ar' = FAr in 
equation (6). Imposing the limits over r in equation (9) recovers the S factor in the numerical integration of 
equation (11) performed in Section 2. Furthermore noting that = ro//3Fc, and defining = (1 + z)e, we 
obtain 

Mt) = (^)-3+" X //[^;ma.(l - ^),min(l - /3m,-, ^)] oc e«i;3+« . (B3) 

za^ '0 '0 '-^z 

The final proportionality holds provided that is in the range satisfying the Heaviside function. When 
a = 0, corresponding to emission at the peak of the vF,, spectrum, f^pki^) oc oc 6^^,. This follows because 
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epk oc t^^, as is apparent by inspecting the limits in the Heaviside function {pctz/ro oc e'^/Tcz, so that 

The validity of eq. (B3) can be checked by deriving the fluence ip = de dt fe{t)/e, using normal- 
ization (B2), from which eq. (A2) is recovered. 



C. Relationship between /^(t) and Uq in the Curvature Limit 

We now derive an approximate analytic expression for a radiation pulse in the curvature limit. Substi- 
tuting expressions (6) and (7) for the comoving spectral energy density into eq. (5), the r-intcgral can be 
approximately solved by letting J drr^[. . -j/Ar' rg|c?r/dr'|[. . .] = Sr^l . .]. In the limit F 1, one obtains 

- ^ 3-a + -Y^^'^ ' (^^^ 

where 

Qa = [ma^(l, ^)]-^+« - [min(4^^«, !^)]-3+a , (C2) 
I + 11 e 

Q, = [max(l, !£M)]-3+;< _ [niin(4r2, , (C3) 

u = 2T^ctz/ro, 7] = TcAt'/ro = ilt/r]r, and we set = —1 (otherwise the terms AF^ are replaced with 
2r^(l — in eqs. (C2) and (C3) above). The vFi, peak energy observed at the start of the pulse is 
denoted by £^^,0 = SrCy^. o/(l + z). 

By examining the limits in eq. (C2), one finds that 



Qa 



( 1 - , i<u<^ <l + r] , 

1 _ (££|^)-3+a ^ 1 < ££to < y < 1 + ^ ^ 

( u \—3+a _ ( <^pk,o \—3+a i ^ u ^ <^pk,o ^ (C4) 



with related expressions for Qi,. In the limit <C 1, corresponding to the curvature limit where variability 
arises principally from curvature effects, the fourth relation in eq. (C4) applies, giving 

/, ^ l£!^2ro£! ^[(^)a„-3+a;j(!£M _ ^) ^ (^)fey-3+6^(e _ »)] . (C5) 

This expression applies equally to the late-time asymptote t ^ (1 + z)At' /2V. At the peak of the uF^ 
spectrum, a = 0, and 

4cM()rgr^ .21^.-3 

recovering the dependence derived in eq. (B3). Note that because Tq oc 7}^, /e oc r?r^t- Eq. (C5) relates f^{t) 



D. Searching for Curvature Effects in GRB Pulses 

If the GRB spectral flux is described by a power-law spectrum with vF^ index a, then curvature effects 
would produce the behavior 

/eW (Di) 

The FWHM width of the curvature spectrum in such a regime is given by = — l]tpk, where 

tpk = (1 + ^)ro/2r^c. Hence the expression tq = 2r'^ctFWHM/{[^^^^^~''^ - + z)} relates the blast 
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wave radius ro to F, given the observables z and tpwHM — provided that the pulse shape is determined by 
curvature effects. Curvature effects also dictate that /e . oc c^^. By examining the variation of intensity as 
a function of Cpk for two GRBs, Soderberg & Fenimore (2001) searched for the signature of shell curvature 
using an expression analagous to eq. (Dl), though without success. In this case, ff oc ^^-{<^+f>) ^ where 
(o + 6) is the mcSiii iiid.6x of th.6 photoii flux wittiiii the intcrvcil contSiining Cpi^ used to measure fe^^ (to first 
order, {a + b) = 0). 

Eq. (C6) can be used to derive an expression relating Cpfc to the photon fluence cp in the curvature limit. 
One simply obtains 

^ = Jl-^, (D2) 

^pk,0 V V'tot 

where Cpkfi refers to the value of epk at the beginning of the pulse, and (ptot refers to the total fluence. This 
expression represents an alternative analytic form to the relation Cpk/ephfl = &^p{—^/^tot) proposed by 
Liang & Kargatis (1996), and also derives from eq. (C5), provided that a and b are independent of time. 

When expression (D2) deviates from observational data, as will often be the case since the approximations 
leading to the curvature pulse are rarely expected to be realized in GRB colliding shells, then curvature effects 
can still be sought by numerically evaluating eqs. (11) or (CI) to obtain more general epk-^ relations. These 
equations can also be used to fit pulse profiles directly. Such an approach would place the phenomenological 
treatments of Koccvski, Ryde, & Liang (2003) and Ryde et al. (2004) on a physical basis, and can be extended 
to treat realistic electron injection and loss processes. Such results can then be compared with predictions 
of the external shock models for the prompt phase. 
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